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Abstract . This report summarizes the research performed at Rice University 
during the period 1984-86 under NASA Grant No. NAG-1-516 on optimal flight 
trajectories in the presence of windshear. With particular reference to the 
take-off problem, the topics covered include: equations of motion, problem 
formulation, algorithms, optimal flight trajectories, advanced guidance schemes, 
simplified guidance schemes, and piloting strategies. 

Key Words . Flight mechanics, take-off, optimal trajectories, optimal control, 
feedback control, windshear problems, sequential gradient-restoration algorithm, 
dual sequential gradient- restoration algorithm, guidance strategies, acceleration 
guidance, gamma guidance, theta guidance, piloting strategies. 
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1 . Introduction 

The objective of this study is the determination of optimal flight 
trajectories in the presence of windshear and the development of guidance 
schemes for near-optimum flight in a windshear. This is a wind characterized 
by sharp change in intensity and direction over a relatively small region of 
space. This problem is important in the take-off and landing of both civilian 
airplanes and military airplanes and is key to aircraft safety. 

The research done under NASA Grant No. NAG-1 -51 6 during the period 1984-86 
is reviewed in Section 2. With particular reference to the take-off problem, 
the topics covered include; equations of motion, problem formulation, algorithms, 
optimal flight trajectories, advanced guidance schemes, simplified guidance 
schemes, and piloting strategies. 

The publications completed or in progress are listed in Section 3. The 
abstracts of the publications completed are given in Section 4. Finally, a 
bibliography is presented in Section 5. 
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2. Research Review 

2.1. Background . During the month of November 1983, Captain W. W. Melvin, 
Delta Airlines and ALFA, approached Dr. A. Miele, Professor of Astronautics 
and Mathematical Sciences at Rice University, asking him to become interested 
in the windshear problem. In a meeting which took place on the campus of Rice 
University, Captain Melvin stated his feeling that considerable research had 
been done on the meteorological, aerodynamic, instrumentation, and stability 
aspects of the windshear problem; however, relatively little had been done on 
the flight mechanics aspects; he felt that a fundamental study was needed in 
order to better understand the dynamic behavior of an aircraft in a windshear. 

In the ensuing discussion. Dr. Miele stated his feeling that the determination 
of good strategies for coping with windshear situations was essentially an 
optimal control problem; that the methods of optimal control theory were needed; 
and that, only after having found optimal control solutions, one could properly 
address the guidance problem. 

As a result of this meeting. Dr. Miele prepared a research proposal 
on the optimization and guidance of flight trajectories in a windshear. The 
proposal was funded in August 1984 by NASA-Langley Research Center, with 
Dr. Miele acting as Principal Investigator and Captain Melvin acting as 
Consultant. Dr. R. L. Bowles of NASA-LRC is Project Monitor. Additional 
funds were subsequently obtained through the sponsorship of Boeing Commercial 
Aircraft Company (BCAC) and Air Line Pilots Association (ALPA). 

During the past two years, the Aero-Astronautics Group of Rice University, 
under the direction of Dr. Miele, has done research on both the optimization 
and the guidance of flight trajectories in a windshear. In the first year, 
the Aero-Astronautics Group studied the optimization aspects of flight 
trajectories (Refs. 1-7 and 14); in the second year, the Aero-Astronautics 
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Group studied the guidance aspects (Refs. 8-13 and 15-19). The main results 
are summarized below. 

2.2. Equations of Motion . In Ref. 1, the equations of motion under 
windshear conditions are derived employing three different coordinate systems: 
the Earth-fixed system; the relative wind-axes system; and the absolute 
wind-axes system. The following assumptions are employed: the aircraft is a 
particle of constant mass; flight takes place in a vertical plane; Newton's 
law is valid in an Earth-fixed system; and the wind flow field is steady. 

For the optimization of flight trajectories, any of the previous 
coordinate systems can be used. However, the relative wind-axes system is to 
be preferred, because the windshear terms appear explicity in the dynamical 
equations. Therefore, use of the relative wind-axes system allows an easier 
physical understanding and interpretation of windshear phenomena. 

2.3. Problem Formulation . In Ref. 2, we employ the equations of 
motion written for the relative wind-axes system. First, we supply an analytical 
description of the forces acting on the aircraft. Next, we supply a description 
of the wind flow field. Generally speaking, the wind flow field is two-dimensional 
however, useful one-dimensionaV models can be developed if one refers to the 
near-the-ground behavior of a microburst. 

With reference to take-off, we assume that the power setting is held at 
the maximum value. Indeed, it is logical to think that, if a plane takes off 
under less-than-ideal weather conditions, a prudent pilot employs the maximum 
thrust. With the power setting held at the maximum value, it is clear that the 
only control is the angle of attack a. To obtain realistic trajectories, 
inequality constraints must be imposed on both a and a. Specifically, the angle 
of attack a is subject to the inequality a _< a*, where a* is a prescribed upper 
bound. In addition, its time derivative a is subject to the double inequality 
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-C £ a _< +C, where C is a prescribed constant. 

Concerning the initial conditions, we assume that the initial state 
is given. Concerning the final conditions, four cases are considered; hence, 
four types of boundary conditions are considered. 

(BCO) The state is free at the final point. 

(BCl) The final value of the path inclination is the same as the 
initial value. 

(BC2) The final values of the velocity and the path inclination are the 
same as the initial values. 

(BC3) The final values of the velocity, the path inclination, and the 
angle of attack are the same as the initial values. Therefore, this case implies 

that, if the initial values correspond to quasi-steady flight, then the final 

1 ) 

values also correspond to quasi-steady flight. i 

Concerning the performance indexes, we consider eight fundamental 
optimization problems. 

(PI). This is a least-square problem involving Ah = h - h^, the difference 
between the flight altitude and a reference altitude, assumed to be a linear 
function of the horizontal distance. 

(P2) This is a least-square problem involving Ay = y - Yp» the 

difference between the relative path inclination and a reference value, assumed 

constant. 

(P3) This is a least-square problem involving Ay^ = - Y d» the 

difference between the absolute path inclination and a reference value, assumed 
constant. 

(P4) This is a minimax problem involving Ah = h - hj^, the difference 
between the flight altitude and a reference altitude, assumed constant. 

(P5) This is a minimax problem involving GAh = G(h - h^), the weighted 

difference between the flight altitude and a reference altitude, assumed constant. 
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as in (P4). Here, G(t) is a prescribed weighting function. 

(P6) This is a minitnax problem involving Ah = h - h^, the 
difference between the flight altitude and a reference altitude, assumed to be 
a linear function of the horizontal distance, as in (PI). 

(P7) This is a minimax problem involving Ay = y - Y|^» tbe difference 
between the relative path inclination and a reference value, assumed constant, 
as in (P2). 

(P8). This is a minimax problem involving Ay„ = " Y„n» the 

difference between the absolute path inclination and a reference value, assumed 
constant, as in (P3). 

Problems (P1)-(P3) are least-square problems of the Bolza type. Problems 
(P4)-(P8) are minimax problems of the Chebyshev type, which can be converted 
into Bolza problems through suitable transformations. Hence, Problems (P1)-(P8) 
are particular cases of the following general problem: 

(P) Minimize a functional with respect to the state vector x(t), the 
control vector u(t), and the parameter vector u which satisfy a system of 
differential constraints, initial constraints, and final constraints. 

2.4. Algorithms . From the previous section, it is clear that one is faced 
with a wide variety of problems of optimal control, depending on the particular 
performance index chosen and the particular type of boundary conditions chosen. 
These problems are further complicated by the presence of inequality constraints 
on the control (a) and the time derivative of the control (a). Therefore, a 
powerful algorithm is necessary to solve the problems under consideration. 

In Ref. 3, we present the algorithm useful for solving Problem (P) 
on a digital computer, more specifically, the sequential gradient-restoration 
algorithm (SGRA). Both the primal formulation and the dual formulation are 
presented. Depending on whether the primal formulation is used or the dual 
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formulation is used, one obtains a primal sequential gradient- restoration 
algorithm (PSGRA) or a dual sequential gradient-restoration algorithm (DS6RA). 

The systems of Lagrange multipliers associated with the gradient phase 
of SGRA and the restoration phase of SGRA are examined. For each phase, 
it is shown that the Lagrange multipliers are endowed with a duality property: 
they minimize a special functional, quadratic in the multipliers, subject to 
the multiplier differential equations and boundary conditions, for given state, 
control, and parameter. These duality properties have considerable computational 
implications: they allow one to reduce the auxiliary optimal control problems 
associated with the gradient phase and the restoration phase of SGRA to 
mathematical programming problems involving a finite number of parameters as 
unknowns. 

Numerical experimentation has shown that, for nonstiff problems of flight 
mechanics, DSGRA is computationally more efficient than PSGRA. In particular, 
for the problems under consideration, DSGRA is about 10% more efficient that 
PSGRA. Hence, the subsequent numerical experiments are based on the use of 
DSGRA. 

2.5. Optimal Trajectories . Optimal trajectories were computed for the 
Boeing B-727 aircraft, using the sequential gradient-restoration algorithm and 
the NAS-AS-9000 computer of Rice University. Among the performance indexes 
(P1)-(P8), the most reliable one was found to be (P8), based on the deviation 
of the absolute path inclination from a reference value. After computing several 
hundred optimal trajectories, certain general conclusions became apparent (see 
Refs. 1-7 and 14): 

(i) the optimal trajectories achieve minimum velocity near the time 


when the shear ends; 
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(ii) the optimal trajectories require an initial decrease in the angle 
of attack, followed by a gradual increase; the maximum permissible angle of 
attack a* (stick-shaker angle of attack) is achieved near the time when the shear 
ends; 

(iii) for weak- to -mode rate windshears,the optimal trajectories are 
characterized by a monotonic climb; the average value of the path inclination 
decreases as the intensity of the shear increases; 

(iv) for relatively severe windshears, the optimal trajectories are 
characterized by an initial climb, followed by nearly horizontal flight, 
followed by renewed climbing after the aircraft has passed through the shear 
region; 

(v) weak-to-moderate windshears and relatively severe windshears are 
survivable employing an optimized flight strategy; however, extremely severe 
windshears are not survivable, even employing an optimized flight strategy; 

(vi) in relatively severe windshears, optimal trajectories have a 
much better survival capability than constant angle of attack trajectories 
(for instance, maximum angle of attack trajectories or maximum lift-to-drag 
ratio trajectories); in addition, they have a better survival capability 
than constant pitch trajectories. 

2.6. Advanced Guidance Schemes . The computation of the optimal trajectories 
requires global information of the wind flow field; that is, it requires the 
knowledge of the wind velocity components at every point of the region of space 
in which the aircraft is flying. In practice, this global information is not 
available; even if it were available, there would not be enough computing 
capability onboard and enough time to process it adequately. As a consequence, 
one must think of optimal trajectories as merely benchmark trajectories that it 
is desirable to approach in actual flight. 
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Since global information is not available, the best that one can do is 
to employ local information on the wind flow field, in particular, local 
information on the wind acceleration and the downdraft. Therefore, the 
guidance problem must be addressed in these terms: Assuming that local information 
is available on the wind acceleration, the downdraft, as well as the state 
of the aircraft, we wish to guide an aircraft automatically or semi automatically 
in such a way that the key properties of the optimal trajectories are preserved. 

Based on the idea of preserving the properties of the optimal trajectories, 
four guidance schemes were developed at Rice University: 

(a) acceleration guidance, based on the relative acceleration; 

(b) absolute gamma guidance, based on the absolute path inclination; 

(c) relative gamma guidance, based on the relative path incliantion; 

(d) theta guidance, based on the pitch attitude angle. 

The details of these guidance schemes are omitted and can be found in 
Refs. 8-13 as well as in Refs. 15-17. 

Guidance trajectories were computed for the Boeing B-727 aircraft using 
the above guidance schemes and the NAS-AS-9000 computer of Rice University. It 
was found that all of the above guidance schemes (in particular, the acceleration 
guidance and the gamma guidance) produce trajectories which are quite close to 
the optimum trajectories. In addition, the resulting near-optimum trajectories 
are superior to the trajectories arising from alternative guidance schemes 
(for instance, constant angle of attack trajectories and constant pitch trajectories). 

2.7. Simplified Guidance Schemes and Piloting Strategies . As stated 
above, the previous advanced guidance schemes require local information on the 
wind flow field and the state of the aircraft. While this information will be 
available in future aircraft, it might not be available on current aircraft. 
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For current aircraft, one way to survive a windshear encounter is to 
use a constant pitch attitude technique; this technique has been advocated by 
many experts on the windshear problem. An alternative, simple technique 
is the quick transition to horizontal flight, based on the properties of the 
optimal trajectories. 

The quick transition to horizontal flight requires an initial decrease 
of the angle of attack, so as to decrease the path inclination to nearly 
horizontal. Then, nearly horizontal flight is maintained during the shear 
encounter. 

For relatively severe windshears, the quick horizontal flight transition 
technique yields trajectories which are competitive with those of the advanced 
guidance schemes discussed previously. In addition, for relatively severe 
windshears, the quick horizontal flight transition technique yields trajectories 
which have better survival capabilities than those associated with other guidance 
techniques, such as constant angle of attack or constant pitch. 

Work on the quick horizontal flight transition technique is nearly 
completed and publication of two reports is expected shortly (Refs. 18-19). 
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4. Abstracts of Publications 

4.1. MIELE, A., WANG, T., and MELVIN, W. W., Optimal Flight Trajectories in 
the Presence of Windshear, Part 1, Equations of Motion , Rice University, 
Aero-Astronautics Report No. 191, 1985. 

Abstract . This report is the first of series dealing with the determination 
of optimal flight trajectories in the presence of windshear. This is a wind 
characterized by sharp change in intensity and direction over a relatively small 
region of space. This problem is important in the take-off and landing of 
both civilian airplanes and military airplanes and is key to aircraft safety. 

It is assumed that: the aircraft is a particle of constant mass; flight 
takes place in a vertical plane; Newton's law is valid in an Earth-fixed system; 
and the wind flow field is steady. 

Under the above assumptions, the equations of motion under windshear 
conditions are derived employing three different coordinate systems: the 
Earth-fixed system; the relative wind-axes system; and the absolute wind-axes 
system. Transformations are supplied which allow one to pass from one system 
to another. 

4.2. MIELE, A., WANG, T., and MELVIN, W. W., Optimal Flight Trajectories in 
the Presence of Windshear, Part 2, Problem Formulation, Take-Off , Rice 
University, Aero-Astronautics Report No. 192, 1985. 

Abstract . This report is the second of a series dealing with the 
determination of optimal flight trajectories in the presence of windshear. We 
employ the equations of motion written for the relative wind-axes system. We 
supply an analytical description of the forces acting on the aircraft, as well 
as a description of the wind flow field. Then, with reference to take-off, 
we formulate eight fundamental optimization problems [Problems (P1)-(P8)] 
under the assumptions that the power setting is held at the maximum value and 
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that the airplane is controlled through the angle of attack. 

Problems (P1)-(P3) are least-square problems of the Bolza type. Problems 
(P4)-(P8) are minimax problems of the Chebyshev type, which can be converted 
into Bolza problems through suitable transformations. Hence, (P1)-(P8) can 
be solved employing the family of sequential gradient-restoration algorithms 
(SGRA), developed for optimal control problems of the Bolza type. 

4.3. MIELE, A., WANG, T., and MELVIN, W. W., Optimal Flight Trajectories in 
the Presence of Windshear, Part 3, Algorithms , Rice University, Aero- 
Astronautics Report No. 193, 1985. 

Abstract . This report is the third of a series dealing with the 
determination of optimal flight trajectories in the presence of windshear. We 
consider the following general problem of the Bolza type [Problem (P)]: Minimize 
a functional with respect to the state vector x(t), the control vector u(t), 
and the parameter vector tr which satisfy a system of differential constraints, 
initial constraints, and final constraints. 

We present the algorithms useful for solving Problem (P) on a digital 
computer, more specifically, sequential gradient-restoration algorithms (SGRA). 
Both the primal formulation and the dual formulation are presented. Depending 
on whether the primal formulation is used or the dual formulation is used, one 
obtains a primal sequential gradient-restoration algorithm (PSGRA) or a dual 
sequential gradient-restoration algorithm (DSGRA). 

The system of Lagrange multipliers associated with the gradient phase 
of SGRA and the restoration phase of SGRA is examined. For each phase, it 
is shown that the Lagrange multipliers are endowed with a duality property: 
they minimize a special functional, quadratic in the multipliers, subject to 
the multiplier differential equations and boundary conditions, for given state, 
control, and parameter. These duality properties have considerable computational 
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implications: they allow one to reduce the auxiliary optimal control problems 
associated with the gradient phase and the restoration phase of S6RA to 
mathematical programming problems involving a finite number of parameters as 
unknowns. 

4.4. MIELE, A., WANG, T., and MELVIN, W. W., Optimal Flight Trajectories in 
the Presence of Windshear, Part 4, Numerical Results, Take-Off , Rice 
University, Aero-Astronautics Report No. 194, 1985. 

Abstract . This report is the fourth of a series dealing with 
the determination of optimal flight trajectories in the presence of windshear. 

We obtain numerical results for the take-off problem, employing the dual 
sequential gradient-restoration algorithm (DSGRA). We investigate a large 
number of combinations of performance indexes, boundary conditions, windshear 
models, and windshear intensities. Inequality constraints are imposed on only 
the angle of attack or on both the angle of attack and its time derivative. 

The following conclusions are reached: (i) optimal trajectories are 
considerably superior to constant angle of attack trajectories; (ii) optimal 
trajectories achieve minimum velocity at about the time when the windshear 
ends; (iii) optimal trajectories can be found which transfer an aircraft from a 
quasi-steady condition to a quasi-steady condition through a windshear; (iv) 
among the optimal trajectories investigated, those minimaximizing |Ay| are 
of particular interest, because the altitude distribution exhibits a monotonic 
behavior; this is true for a moderate windshear and a relatively severe windshear 
(v) an extremely severe windshear cannot be survived, even employing an optimized 
flight strategy; and (vi) the sequential gradient-restoration algorithm has 
proven to be a powerful algorithm for solving the problem of the optimal flight 
trajectories in a windshear. 
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4.5. MIELE, A., Summary Report on NASA Grant No. NAG-1 -51 6, Optimal Flight 
Trajectories in the Presence of Windshear, 1985-86 , Rice University, 
Aero-Astronautics Report No. 195, 1985. 

Abstract . This report summarizes the research performed at Rice University 
during the period 1984-85 under NASA Grant No. NAG-1 -51 6 on optimal flight 
trajectories in the presence of windshear. The topics covered include: equations 
of motion, problem formulation (take-off), algorithms, and numerical results 
(take-off) . 

4.6. MIELE, A., WANG, T., and MELVIN, W. W.. Optimal Take-Off Trajectories in 
the Presence of Windshear , Paper No. AIAA-85-1843-CP, AIAA Atmospheric 
Flight Mechanics Conference, Snowmass, Colorado, 1985. 

Abstract . This paper is concerned with optimal flight trajectories in 
the presence of windshear. With particular reference to take-off, eight fundamental 
optimization problems [Problems (P1)-(P8)] are formulated under the assumptions 
that the power setting is held at the maximum value and that the airplane is 
controlled through the angle of attack. 

Problems (P1)-(P3) are least-square problems of the Bolza type. Problems 
(P4)-(P8) are minimax problems of the Chebyshev type, which can be converted 
into Bolza problems through suitable transformations. These problems are 
solved employing the dual sequential gradient-restoration algorithm (DSGRA) for 
optimal control problems. 

Numerical results are obtained for a large number of combinations of 
performance indexes, boundary conditions, vnndshear models, and windshear 
intensities. However, for the sake of brevity, the presentation of this paper 
is restricted to Problem (P6), minimax |Ah|, and Problem (P7), minimax [Ay]. 
Inequality constraints are imposed on the angle of attack and the time derivative 
of the angle of attack. 
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The following conclusions are reached: (i) optimal trajectories are 
considerably superior to constant angle of attack trajectories; (ii) optimal 
trajectories achieve minimum velocity at about the time when the windshear ends; 
(iii) optimal trajectories can be found which transfer an aircraft from a 
quasi-steady condition to a quasi-steady condition through a windshear; (iv) 
among the optimal trajectories investigated, those minimaximizing |Ay| are of 
particular interest, because the altitude distribution exhibits a monotonic 
behavior; this is true for a moderate windshear and a relatively severe windshear 
(v) an extremely severe windshear cannot be survived, even employing an optimized 
flight strategy; and (vi) the sequential gradient-restoration algorithm (SGRA) 
has proven to be a powerful algorithm for solving the problem of the optimal 
flight trajectories in a windshear. 

4.7. MIELE, A., WANG, T., and MELVIN, W. W., Optimal Take-Off Trajectories 
in the Presence of Windshear , Journal of Optimization Theory and 
Applications, Vol. 49, Mo. 1, pp. 1-45, 1986. 

Abstract . This paper is concerned with optimal flight trajectories in the 
presence of windshear. With particular reference to take-off, eight fundamental 
optimization problems [Problems (P1)-(P8)] are formulated under the assumptions 
that the power setting is held at the maximum value and that the airplane is 
controlled through the angle of attack. 

Problems (P1)-(P3) are least-square problems of the Bolza type. Problems 
(P4)-(P8) are minimax problems of the Chebyshev type, which can be converted 
into Bolza problems through suitable transformations. These problems are 
solved employing the dual sequential gradient-restoration algorithm (DSGRA) for 
optimal control problems. 

Numerical results are obtained for a large number of combinations of 
performance indexes, boundary conditions, windshear models, and windshear 
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intensities. However, for the sake of brevity, the presentation of this paper 
is restricted to Problem (P6), minimax lAh|, and Problem (P7), minimax |Ay|. 
Inequality constraints are imposed on the angle of attack and the time derivative 
of the angle of attack. 

The following conclusions are reached: (i) optimal trajectories are 
considerably superior to constant-angle-of -attack trajectories; (ii) optimal 
trajectories achieve minimum velocity at about the time when the windshear ends; 
(iii) optimal trajectories can be found which transfer an aircraft from a 
quasi-steady condition to a quasi-steady condition through a windshear; (iv) 
as the boundary conditions are relaxed, a higher final altitude can be achieved, 
albeit at the expense of a considerable velocity loss; (v) among the optimal 
trajectories investigated, those solving Problem (P7) are to be preferred, 
because the altitude distribution exhibits a monotonic behavior; in addition, 
for boundary conditions BC2 and BC3, the peak angle of attack is below the 
maximum permissible value; (vi) moderate windshears and relatively severe 
windshears are survivable employing an optimized flight strategy; however, 
extremely severe windshears are not survivable, even employing an optimized 
flight strategy; and (vii) the sequential gradient-restoration algorithm (SGRA), 
employed in its dual form (DSGRA), has proven to be a powerful algorithm for 
solving the problem of the optimal flight trajectories in a windshear, 

4.8. MIELE, A., WANG, T., and MELVIN, W. W., Guidance Strategies for Near- 

Optimum Performance in a Windshear, Part 1 , Take-Off, Basic Strategies , 

Rice University, Aero-Astronautics Report No. 201, 1986. 

Abstract . This paper is concerned with guidance strategies for near- 
optimum performance in a windshear. The take-off problem is considered with 
reference to flight in a vertical plane. 



19 


AAR-203 


First, trajectories for optimum performance in a windshear are determined 
for different windshear models and different windshear intensities. Use is 
made of the methods of optimal control theory in conjunction with the dual 
sequential gradient-restoration algorithm (DSGRA) for optimal control problems. 

In this approach, global information on the wind flow field is needed. 

Then, guidance strategies for near-optimum performance in a windshear 
are developed, starting from the optimal trajectories. Specifically, three 
guidance schemes are presented; (A) gamma guidance, based on the relative 
path inclination; (B) theta guidance, based on the pitch attitude angle; and 
(C) acceleration guidance, based on the relative acceleration. In this approach, 
local information on the wind flow field is needed. 

Numerical experiments show that guidance schemes (A), (B), (C) produce 
trajectories which are quite close to the optimum trajectories. In addition, 
the near-optimum trajectories are considerably superior to the trajectories 
arising from alternative guidance schemes. 

An important characteristic of guidance schemes (A), (B), (C) is their 
simplicity. Indeed, these guidance schemes are implementable using available 
instrumentation and/or modification of available instrumentation. 

4.9. MIELE, A., WANG, T., and MELVIN, W. W., Guidance Strategies for Near- 

Optimum Performance in a Windshear, Part 2, Take-Off, Comparison Strategies , 
Rice University, Aero-Astronautics Report No. 202, 1986. 

Abstract . In the previous paper, trajectories for optimum performance 
in a windshear were determined for different windshear models and different 
windshear intensities. Then, guidance strategies for near-optimum performance 
in a windshear were developed, starting from the optimal trajectories. Specifically, 
three guidance schemes were presented: (A) gamma guidance, based on the relative 
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path inclination; (B) theta guidance, based on the pitch attitude angle; and 
(C) acceleration guidance, based on the relcltive acceleration. 

In this report, several comparison strategies are investigated for the 
sake of completeness, more specifically: (D) constant alpha guidance; (E) constant 
velocity guidance; (F) constant theta guidance; (G) constant relative path 
inclination guidance; (H) constant absolute path inclination guidance; and 
(I) linear altitude distribution guidance. 

Numerical experiments show that guidance schemes (D) through (I) are 
inferior to guidance schemes (A) through (C) for a variety of technical reasons. 

In particular, for the case where the wind velocity difference is AW = 120 ft sec 

X 

and the windshear intensity is AW /Ax = 0.030 sec“\ it is shown that the 
trajectories associated with guidance schemes (D) through (I) hit the ground, 
while the trajectories associated with guidance schemes (A) through (C) clear 
the ground. 

4.10. MIELE, A., WANG, T., and MELVIN, W. W., Guidance Strategies for Near- 
Optimum Take-Off Performance in a Windshear , Paper No. AIAA-86-0181 , 

AIAA 24th Aerospace Sciences Meeting, Reno, Nevada, 1986. 

Abstract . This paper is concerned with guidance strategies for near- 
optimum performance in a windshear. The take-off problem is considered with 
reference to flight in a vertical plane. 

First, trajectories for optimum performance in a windshear are 
determined for different windshear models and difference windshear intensities. 

Use is made of the methods of optimal control theory in conjunction with the 
dual sequential gradient-restoration algorithm (DSGRA) for optimal control 
problems. In this approach, global information on the wind flow field is needed. 
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Then, guidance strategies for near-optimum performance in a windshear 
are developed, starting from the optimal trajectories. Specifically, three 
guidance schemes are presented: (A) gamma guidance, based on the relative 
path inclination; (B) theta guidance, based on the pitch attitude angle; and 
(C) acceleration guidance, based on the relative acceleration. In this approach, 
local information on the wind flow field is needed. 

Numerical experiments show that guidance schemes (A), (B), (C) produce 
trajectories which are quite close to the optimum trajectories. In addition, 
the near-optimum trajectories are considerably superior to the trajectories 
arising from alternative guidance schemes. 

An important characteristic of guidance schemes (A), (B), (C) is their 
simplicity. Indeed, these guidance schemes are implementable using available 
instrumentation and/or modification of available instrumentation. 

4.11. MIELE, A., WANG, T., and MELVIN, W. W.. Guidance Strategies for Near- 
Optimum Take-Off Performance in a Windshear , Journal of Optimization 
Theory and Applications, Vol. 50, No. 1, pp. 1-47, 1986. 

Abstract . This paper is concerned with guidance strategies for near- 
optimum performance in a windshear. The take-off problem is considered with 
reference to flight in a vertical plane. 

First, trajectories for optimum performance in a windshear are 
determined for different windshear models and different windshear intensities. 

Use is made of the methods of optimal control theory in conjunction with the 
dual sequential gradient-restoration algorithm (DSGRA) for optimal control 
problems. In this approach, global information on the wind flow field is needed. 

Then, guidance strategies for near-optimum performance in a windshear 
are developed, starting from the optimal trajectories. Specifically, three 
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guidance schemes are presented: (A) ganma guidance, based on the relative 
path inclination; (B) theta guidance, based on the pitch attitude angle; and 
(C) acceleration guidance, based on the relative acceleration. In this 
approach, local information on the wind flow field is needed. 

Next, several alternative schemes are investigated for the sake of 
completeness, more specifically: (D) constant alpha guidance; (E) constant 
velocity guidance; (F) constant theta guidance; (G) constant relative path 
inclination guidance; (H) constant absolute path inclination guidance; and 
(I) linear altitude distribution guidance. 

Numerical experiments show that guidance schemes (A)-(C) produce 
trajectories which are quite close to the optimum trajectories. In addition, 
the near-optimum trajectories associated with guidance schemes (A)-(C) are 
considerably superior to the trajectories arising from the alternative guidance 
schemes (D)-(I). 

An important characteristic of guidance schemes (A)-(C) is their 
simplicity. Indeed, these guidance schemes are implementable using available 
instrumentation and/or modification of available instrumentation. 

4.12. MIELE, A., WANG, T., and MELVIN, W. W., Optimization and Acceleration 

Gui dance of FI i qht T ra jectOri es i h a Wi ndshear , Paper No. AIAA-86-2036-CP, 
AIAA Guidance, Navigation, and Control Conference, Williamsburg, Virginia, 
1986. 

Abstract. This paper is concerned with guidance strategies for near- 
optimum performance in a windshear. The take-off problem is considered with 
reference to flight in a vertical plane. In addition to the horizontal shear, 
the presence of a downdraft is assumed. 

First, trajectories for optimum performance in a windshear are determined 
for different windshear models and different windshear intensities. Use is 
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made of the methods of optimal control theory in conjunction with the dual 
sequential gradient-restoration algorithm (DSGRA) for optimal control problems. 

In this approach, global information on the wind flow field is needed. 

Then, guidance strategies for near-optimum performance in a windshear 
are developed, starting from the optimal trajectories. Specifically, an 
acceleration guidance scheme, based on the relative acceleration, is presented 
in both analytical form and feedback control form. In this approach, local 
information on the wind flow field is needed. 

Numerical experiments show that the acceleration guidance scheme produces 
trajectories which are quite close to the optimum trajectories. In addition, 
the near-optimum trajectories are considerably superior to the trajectories 
arising from alternative guidance schemes. 

An important characteristic of the acceleration guidance scheme is its 
simplicity. Indeed, this guidance scheme is implementable using available 
instrumentation and/or modification of available instrumentation. 

4.13. MIELE, A., WANG, T., and MELVIN. W. W., Optimization and Gamma/Theta 
Guidance of Flight Trajectories in a Windshear , Paper No. ICAS-86-564, 
15th Congress of the International Council of the Aeronautical Sciences, 
London, England, 1986. 

Abstract . This paper is concerned with guidance strategies for near- 
optimum performance in a windshear. The take-off problem is considered with 
reference to flight in a vertical plane. In addition to the horizontal shear, 
the presence of a downdraft is assumed. 

First, trajectories for optimum performance in a windshear are determined 
for different windshear models and different windshear intensities. Use is 
made of the methods of optimal control theory in conjunction with the dual 
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sequential gradient-restoration algorithm (DSGRA) for optimal control problems. 

In this approach, global information on the wind flow field is needed. 

Then, guidance strategies for near-optimum performance in a windshear 
are developed, starting from the optimal trajectories. Specifically, three 
guidance schemes are presented: the absolute ganma guidance scheme, based 
on the absolute path inclination; the relative gamma guidance scheme, based on 
the relative path inclination; and the theta guidance scheme, based on the 
pitch attitude angle. In this approach, local information on the wind flow 
field is needed. 

Numerical experiments show that the gairana/ theta guidance schemes produce 
trajectories which are quite close to the optimum trajectories. In addition, 
the near-optimum trajectories are considerably superior to the trajectories 
arising from alternative guidance schemes. 

An important characteristic of the gamma/ theta guidance schemes is their 
simplicity. Indeed, these guidance schemes are implementable using available 
instrumentation and/or modification of available instrumentation. 

4.14. MIELE, A., WANG, T., and MELVIN, W. W., Quasi -Steady Flight to Quasi- 
Steady Flight Transition in a Windshear: Trajectory Optimization , 6th 
I FAC Workshop on Control Applications of Nonlinear Programming and 
Optimization, London, England, 1986. 

Abstract . This paper is concerned with the optimal transition of an 
aircraft from quasi -steady flight to quasi -steady flight in a windshear. The take-off 
problem is considered with reference to flight in a vertical plane. In addition 
to the horizontal shear, the presence of a downdraft is considered. It is 
assumed that the power setting is held at the maximum value and that the 
aircraft is controlled through the angle of attack. Inequality constraints 
are imposed on both the angle of attack and its time derivative. 
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The optimal transition problem is formulated as a minimax problem or 
Chebyshev problem of optimal control: the performance index being minimized 
is the peak value of the modulus of the difference between the absolute path 
inclination and a reference value, assumed constant. By suitable transformation, 
the Chebyshev problem is converted into a Bolza problem. Then, the Bolza 
problem is solved employing the dual sequential gradient- restoration algorithm 
(DSGRA) for optimal control problems. 

Numerical experiments are performed for different windshear intensities 
and different windshear models. Three basic windshear models are considered: 
(WSl) it includes the horizontal shear and neglects the downdraft; (WS2) it 
neglects the horizontal shear and includes the downdraft; (WS3) it includes 
both the horizontal shear and the downdraft. 

The numerical results lead to the following conclusions: (i) not only the 
transition from quasi -steady flight to quasi-steady flight in a windshear is 
possible, but it can be performed in an optimal way; (ii) for weak-to-moderate 
shear/downdraft combinations, the optimal transition is characterized by a 
monotonic climb; in this monotonic climb, the absolute path inclination is 
nearly constant through the shear region; this constant value decreases as 
the shear/downdraft intensity increases; (iii) for severe shear/downdraft 
combinations, the optimal transition is characterized by an initial climb, 
followed by nearly horizontal flight in the shear region, followed by 
renewed climbing in the after-shear region; (iv) the relative velocity decreases 
in the shear region, achieves a minimum value at about the end of the shear, 
and then increases in the after-shear region; (v) in the after-shear region, 
the absolute path inclination increases at a rapid rate after the velocity 
recovery is almost completed; and (vi) the dual sequential gradient-restoration 
algorithm has proved to be a powerful algorithm for solving the problem of 
the optimal transition in a windshear. 


26 


AAR-203 


5. Bibliography 

Flight Mechanics, Applied Aerodynamics, Stability and Control 

5.1. ABZUG, M. J., Airspeed Stability under Windshear Conditions , Journal 
of Aircraft, Vol. 14, 1977. 

5.2. FROST, W., and CROSBY, B., Investigations of Simulated Aircraft 
Flight through Thunderstorm Outflows , NASA, Contractor Report 
No. 3052, 1978. 

5.3. FROST, W., and RAVIKUMAR, R., Investigation of Aircraft Landing 
in Variable Wind Fields , NASA, Contractor Report No. 3073, 1978. 

5.4. FROST, W., Flight in Low-Level Windshear , NASA, Contractor Report 
No. 3678, 1982. 

5.5. GERA, J., The Influence of Vertical Wind Gradients on the 
Longitudinal Motion of Airplanes , NASA, Technical Note No. D-6430, 
1971. 

5.6. GERA, J., Longitudinal Stability and Control in Windshear with 
Energy Height Rate Feedback , NASA, Technical Memorandum No. 81828, 
1980. 

5.7. HAINES, P. A., and LUERS, J. K., Aerodynamic Penalties of Heavy 

Rain on a Landing Aircraft , NASA, Contractor Report No. 156885, 1982. 

5.8. HOUBOLT, J., Survey on Effect of Surface Winds on Aircraft Design 
and Operations and Recommendations for Needed Wind Research , NASA, 
Contractor Report No. 2360, 1973. 

5.9. MCCARTHY, J., BLICK, E. F., and BENSCH, R. R., Jet Transport 
Performance in Thunderstorm Windshear Conditions, NASA, Contractor 


Report No. 3207, 1979. 


27 


AAR-203 


5.10. MCCARTHY, J., and NORVIEL, V., Numerical and Flight Simulator 
Test of the Flight Deterioration Concept , NASA, Contractor Report 
No. 3500, 1982. 

5.11. MELVIN, W. W., Effects of Windshear on Approach with Associated 
Faults of Approach Couplers and Flight Directors , Paper Presented 
at the AIAA Aircraft Design and Operations Meeting, Los Angeles, 
California, 1969. 

5.12. MELVIN, W. W., The Dynamic Effect of Windshear . Pilot Safety Exchange 
Bulletin, Flight Safety Foundation, Arlington, Virginia, 1975. 

5.13. MELVIN, W. W., What You Don't Know about Windshear Can Kill You , 

Flight Crew, Vol. 2, 1980. 

5.14. REEVES, P. M., et al. Development and Applications of a Nongaussian 
Atmospheric Turbulence Model for Use in Flight Simulators , NASA, 
Contractor Report No. 2451, 1974. 

5.15. TURKEL, B. S., and FROST, W., Pilot-Aircraft System Response to 
Windshear , NASA, Contractor Report No. 3342, 1980. 

5.16. TOBIASON, A. R., NASA Overview on Windshear Research , Oral Presentation, 
National Research Council, Washington, DC, 1983. 

5.17. WANG, S. T., and FROST, W., Atmospheric Turbulence Simulation Techniques 
with Application to Flight Analysis , NASA, Contractor Report No. 3309, 
1980. 

Meteorology 

5.18. CARACENA, F., FLUECK, J. A., and MCCARTHY, 0., Forecasting the 
Likelihood of Microbursts along the Front Range of Colorado , Paper 
Presented at the 13th AMS Conference on Severe Local Storms, Tulsa, 
Oklahoma, 1983. 


28 


AAR- 203 


5.19. DANIEL, J., and FROST, W., 1981 Current Research on Aviation 
Weather (Bi b1 iography) , NASA, Contractor Report No. 3527, 1982. 

5.20. FUJITA, T. T., Downbursts and Microbursts: An Aviation Hazard , 

Paper Presented at the 9th AMS Conference on Radar Meteorology, 

Miami, Florida, 1980. 

5.21. FUJITA, T. T., and CARACENA, F., An Analysis of Three Weather- 
Related Aircraft Accidents , Bulletin of the American Meteorological 
Society, Vol. 58, 1977. 

5.22. GOFF, R. C., Vertical Structure of Thunderstorm Outflows , Monthly 
Weather Review, Vol. 104, 1976. 

5.23. GOFF, R. C., Some Observations of Thunderstorm- Induced Low-Level 
Wind Variations , Journal of Aircraft, Vol. 14, 1977. 

5.24. LEHMAN, J. M., HEFFLEY, R. K., and CLEMENT, W. F., Simulation 

and Analysis of Windshear Hazard , Federal Aviation Administration, 
Report No. FAA-RD-78-7, 1977. 

5.25. PIELKE, R. A., Mesoscale Dynamic Modeling , Advances in Geophysics, 
Academic Press, New York, New York, Vol. 23, 1981. 

5.26. PIELKE, R. A., The Role of Mesoscale Numerical Models in Very-Short- 
Range Forecasting , Nowcasting, Academic Press, New York, New York, 

Vol. XX, 1982. 

5.27. SOWA, D. F., Low-Level Windshear , DC Flight Approach, McDonnell -Douglas 
Aircraft Company, Long Beach, California, No. 20, 1974. 

5.28. SOWA, D. F., The Effect of Terrain Near Airports on Significant 
Low-Level Windshear , Paper Presented at the AIAA Aircraft Systems 
and Technology Conference, Seattle, Washington, 1977. 

5.29. TOWNSEND, J. W., Jr., et al , Low-Altitude Windshear and Its Hazard 
to Aviation, National Academy Press, Washington, DC, 1983. 


29 


AAR-203 


5.30. WURTELE, M. G., Meteorological Conditions Surrounding the Paradise 
Airline Crash of March 1964 , Journal of Applied Meteorology, 

Vol . XX, 1970. 

Optimal Control 

5.31. PONTRYAGIN, L. S., BOLTYANSKII, V. G., GAMKRELIDZE, R. V., and 
MISHCHENKO, E. F., The Mathematical Theory of Optimal Processes , 

John Wiley and Sons (Interscience Publishers), New York, New York, 1962. 

5.32. HESTENES, M. R., Calculus of Variations and Optimal Control Theory , 

John Wiley and Sons, New York, New York, 1966. 

5.33. BRYSON, A. E., Jr., and HO, Y. C., Applied Optimal Control , Blaisdell 
Publishing Company, Waltham, Massachusetts, 1969. 

5.34. LEITMANN, G., The Calculus of Variations and Optimal Control , Plenum 
Publishing Corporation, New York, New York, 1981. 

Numerical Algorithms 

5.35. MIELE, A., PRITCHARD, R. E., and DAMOULAKIS, J. N., Sequential Gradient- 
Restoration Algorithm for Optimal Control Problems , Journal of 
Optimization Theory and Applications, Vol. 5, 1970. 

5.36. MIELE, A., DAMOULAKIS, J. N., CLOUTIER, J. R., and TIETZE, J. L., 
Sequential Gradient-Restoration Algorithm for Optimal Control Problems 
with Nondifferential Constraints , Journal of Optimization Theory and 
Applications, Vol. 13, 1974. 

5.37. MIELE, A., Recent Advances in Gradient Algorithms for Optimal Control 
Problems , Journal of Optimization Theory and Applications, Vol. 17, 

1975. 

5.38. GONZALEZ, S., and MIELE, A., Sequential Gradient-Restoration Algorithm 
for Optimal Control Problems with General Boundary Conditions , Journal 
of Optimization Theory and Applications, Vol. 26, 1978. 


30 


AAR-203 


5.39. MIELE, A., Gradient Algorithms for the Optimization of Dynamic 
Systems , Control and Dynamic Systems, Advances in Theory and 
Application, Edited by C. T. Leondes, Academic Press, New York, 

New York, Vol. 16, 1980. 

5.40. MIELE, A., KUO, Y. M., and COKER, E. M., Modified Quasilinearization 
Algorithm for Optimal Control Problems with Nondifferential Constraints 
and General Boundary Conditions, Parts 1 and 2 , Rice University, 
Aero-Astronautics Reports Nos. 161 and 162, 1982. 

5.41. MIELE, A., BASAPUR. V. K., and COKER, E. M., Combined Gradient- 
Restoration Algorithm for Optimal Control Problans with Nondifferential 
Constraints and General Boundary Conditions, Parts 1 and 2 , Rice 
University, Aero-Astronautics Reports Nos. 163 and 164, 1983. 

Minimax Problems 

5.42. WARGA, J., Minimax Problems and Unilateral Curves in the Calculus 
of Variations , SIAM Journal on Control, Vol. 3, 1965. 

5.43. JOHNSON, C. D., Optimal Control with Chebyshev Minimax Performance 
Index , Journal of Basic Engineering, Vol. 89, 1967. 

5.44. MICHAEL, G. J., Computation of Chebyshev Optimal Control , AIAA Journal, 
Vol. 9, 1971. 

5.45. POWERS, W. F., A Chebyshev Minimax Technigue Oriented to Aerospace 
Trajectory Optimization Problems , AIAA Journal, Vol. 10, 1972. 

5.46. HOLMAKER, K., A Minimax Optimal Control Problem , Journal of Optimization 
Theory and Applications, Vol. 28, 1979. 

5.47. H0Lf*lAKER, K., A Property of an Autonomous Minimax Optimal Control Problem , 
Journal of Optimization Theory and Applications, Vol. 32, 1980. 



31 


AAR-203 


5.48. MIELE, A., MOHANTY, B. P., VENKATARAMAN, P., and KUO, Y. M., Numerical 
Solution of Minimax Problems of Optimal Control, Parts 1 and 2 , 

Journal of Optimization Theory and Applications, Vol . 38, 1982. 
Transformation Techniques 

5.49. VALENTINE, F. A., The Problem of Lagrange with Differential Inequalities 
as Added Side Conditions , Contributions to the Calculus of Variations, 
1933-1937, The University of Chicago Press, Chicago, Illinois, 1937. 

5.50. JACOBSON, D. H., and LELE, M. M., A Transformation Technique for Optimal 
Control Problems with a State Variable Inequality Constraint , IEEE 
Transactions on Automatic Control, Vol. AC-14, 1969. 

5.51. PARK, S. K., A Transformation Method for Constrained Function Minimization , 
NASA, Technical Note No. D-7983, 1975. 

5.52. MIELE, A., WU, A. K., and LIU, C. T., A Transformation Technique for 
Optimal Control Problems with Partially Linear State Inequality 
Constraints , Journal of Optimization Theory and Applications, Vol. 28, 

1979. 

5.53. MIELE, A., and VENKATARAfiAN, P., Optimal Trajectories for Aeroassisted 
Orbital Transfer , Acta Astronautica, Vol. 11, 1984. 

5.54. MIELE, A., and VENKATARAMAN, P., Minimax Optimal Control and Its 
Application to the Reentry of a Space Glider , Recent Advances in the 
Aerospace Sciences, Edited by C. Casci, Plenum Publishing Corporation, 

New York, New York, 1985. 

5.55. MIELE, A., BASAPUR, V. K., and MEASE, K. D., Nearly-Grazing Optimal 
Trajectories for Aeroassisted Orbital Transfer , Journal of the 


Astronautical Sciences, Vol. 33, 1985. 


